Analysis on the imaging properties of a left-handed material slab 
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We investigate in tliis paper the imaging properties of an absorptive left-lianded material (LHM) 
slab. For a line source, a geometric explanation to the reason of the thickness limitation on an ideal 
lossless slab is given. For a lossy slab, the imaging properties are determined by the wavelength, the 
slab thickness, the distance from the source to the nearer boundary of the slab, and the absorption 
Oh' effect. Varying the ratios between these quantities, the image width can be changed from wavelength 

to subwavelength scale. In the former situation, the energy density is mainly concentrated at the 
two image spots. In the later case, though image of subwavelength width appears on the focal plane, 
QQ ■ however, most energy is located at about the two boundaries of the slab. The relations between the 

subwavelength imaging and uncertainty principle are also discussed. 



PACS numbers: 78.20.Ci, 42.30.Wb, 73.20.Mf, 78.66.Bz 



I. INTRODUCTION 
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O ■ Negative refraction of electromagnetic waves by a left-handed material (LHM), first proposed in 1960s by Veselago 
'Si Qi has attracted strong research interests j^lE^IEIll and generated heated debatei?! >8.. ,9.. .10„ ,11.]. Among all the 
J>~>. phenomena that could happen in a LHM, the most fascinating one may be the possibility of "superlensing effect" 
' proposed by Pendry 0; that is, a slab made of uniform and isotropic LHM with both the permittivity e = — 1 
and the permeability fi = —1 acquires a negative refractive index n = —1, which makes this slab a perfect lens. It 
can capture both the propagating and the evanescent waves emitted from a point source placed in front of the slab 
OsJ , and refocuses them into two point images, one inside and the other behind the slab. 

J> ■ Recently, this superlensing effect was questioned by a number of authors 0, |M IS ^i 0| • In Ref. , the authors 
augured that negative refraction of energy flow im plie s the violation of causality principle, and a little amount of 

■ absorption will largely deform the waves. In Ref. |l0|. the authors showed that although there is amplification of 
I ' evanescent waves in an ideal lossless left-handed medium, however, to avoid the divergence of the field energy inside 
I the lens, it must be limited to a thickness smaller than the distance between the line source and the nearer boundary 

■ of the slab, thus perfect imaging is impossble. In addition, a little absorption may destroy the negative refraction 
, effect completely. It was then found that to make a left-handed material physically realizable, the medium must be 

■ dispersive or absorptive. In Ref. Tl'l, the recovery rate for a lossy slab was studied, and the author showed that 
O , the image quality can be significantly affected by the absorption effect. In Ref. 0, the authors showed that 

' ^ . the energy flow indeed goes to the "negative way" when passing through the surface of an absorptive and dispersive 
J>-j' LHM. In I a slab lens of photonic crystal was considered, and the simulation showed that negative refraction of 
. energy flow does not contradict the causality principle. Further in [T5l.lT6i|. the concept of "constant frequency curves" 
introduced in were used to study the refraction behavior of the waves in the the medium. Most interestingly, in 
J> , an all-angle negative refraction photonic crystal slab lens was designed to focus the light into a subwavelength 

' region. 

Although the focusing effect of a LHM slab lens has already been studied by a number of authors, however, in 
most previous studies researchers used some Finite-Difference- Time-Domain (FDTD) method. The method is easy 
■ ■ ' to implement but the physical meanings of the simulation results are not easy to be extracted. In some other studies 
the authors used frequency domain method, however, they usually considered only one single Fourier component of 
the fields. To get a definite result, one has to sum over these Fourier components. 

In this paper we study the imaging problem using a spectrum decomposition method. We first decompose the 
cylindrical wave emitted by a line source into a series of plane waves of different transverse wave numbers. By 
considering the boundary conditions at the source point and the two boundaries of the slab lens, we then can 
determine the transmission and reflection coefficients for each plane wave. These quantities are utilized to construct 
the field function in every space region. 

Our method does not adopt complicated numerical skills, thus makes us easier to get the physical insight. We 
also give a very simple geometrical explanation to the reason of the thickness limitation for the ideal slab lens (the 
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e = fj, = n = —1 case ) [13. Finally, we found that the imaging mechanism for a negative refraction lens system is 
subtler than that of the conventional lens system. 

II. MODEL AND METHOD 

We first describe the setup of the slab system. In this paper we consider only the E-polarized wave, which means 
that the wave propagation direction is parallel to the XZ plane. The x axis is parallel to the two boundaries of the 
slab, and the boundary near the source is the z — plane. A current line source J(r)e~"^* = yJod{r — Tq) e~*"* 
located at Tq = (xo,zo) = (0, zo)j zq < 0, emits monochromatic waves of angular frequency w, thus both the E and 

H fields get a time factor e~*"*. The E field wave radiated from it is Ero(i(r)e~'"* = yAo HQ^\k\r - to\) e~*"*, which 
satisfies 

(V^ + k') E™d(r) = -z^J(r). (1) 

Here Hq^^ (x) is the zeroth order Hankel function of the first kind, Jo and ^0 — — ttlj Jo / are two constants propotional 
to each other, r is the observation point, and k — uj/c and c are the wave number of the cylindrical wave and the 
speed of light in vacuum (outside of the slab), respectively. 

To calculate the total E(r) field, we first introduce the Green's function satisfying 



(V2 + e{z)) G(r,r') = -5^'^\v - r'), (2) 



then the E field is given by 



E(r) = j dVG(r,r')J(r') 



= t^-JoG r,ro (3) 

Here k'^{z) = fc^ = uj^/c^ in the regions outside the slab, and fc^(z) = ejiuP/c^ li Q < z < d. e and /i are the 
permmitivity and permeaility in the slab, respectively. 

To proceed further, the waves have to be decomposed into various Fourier components 5j. Each component has a 
definite kx- It is a plane wave with either a real k^ = yjuj"^ /c^ — k'^, if \kx\ < ^/c, or an imaginary k^, if \kx\ > uj/c. 
In the former case we have a propagating wave, and in the later case the wave is evanescent. 

Write G(r, Tq) as 

G(r,ro) = — / dkxe'''-''g{z,kx), (4) 



then we have 



— t: ~h k (z ) k^ 
dz^ 



g{z,kx) ^ -S{z ~ zq), (5) 



which leads to the boundary condition for g at z — zq: 

g'{z, kx)\z„+ - g'{z, kx)\zo- = -1. (6) 

The continunity conditions for the tangential components of the E and H fields at the two boundaries of the slab lead 
to 

g{,Z kx^\outside — g(.Z ^ kx^\inside : ('^) 
g i^Z , kx)\outaide — 5 {.Z , kx)\inside- (8) 

Define 



Ko \Jk? - fc2, K = y/k'^e^ - fc2, (9) 
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the solution for g is given by 



f e'«ol^-'ol-|-fle'''o"^o'~''' z <0 
— 2iKo ' 



9= S 



giKol^ol T [cos K(z~d)+iJ^ sinK{z-d)] 



. 0<z<d (10) 



7^ gi«o(^-ci+l^ol) 



Z 



> d. 



Here 



(11) 



., . , sin K,d 

R= ^^^^7 (12) 

COSKd- ^ I -J^ + sin Kd 



cosKd-^ (-^ + ^) sinAtd 

and 

I i_ ( K I 

are the transmission and reflection coefficients, calculated from the transfer matrix method pjllr^ 



III. AN IDEAL SLAB 

We now turn to the discussion of an ideal slab lens. For an ideal slab we mean that we can find a frequency wo 
such that for a dispersive medium slab lens medium with frequency dependent permmitivity e(w) and permeability 
/i(w) and zero absorption effect we have s{loo) — ^{loq) = —1. Pendry pointed out in [2] that a slab lens of this kind 
is a perfect lens with n — —1. It focuses the propagating waves and amplifies the evanescent waves, thus can recover 
all the information carried by the wave emitted from the line source. Although Pendry in his derivation showed that 
for a single Fourier component the lens indeed amplifies the evanescent wave and thus the amplitude of the wave can 
be completely recovered, however, he did not sum over these Fourier components to get a result of the total field. In 
[10], the authors showed that if the thickness d of the lens is greater than di — — zqi then the total field will diverge 
inside of the lens. On the other hand, if d < di, there will be no image at all. Thus perfect imaging is impossible. 

Although the thickness limitation discussed in |0 for an ideal LHM slab lens is correct, however, it is hard to 
believe that there is some physical principle that can restrict the slab thickness, if a thinner one can be made. To 
resolve this puzzle, here we give a simple geometrical explanation to the reason of this restriction (See Fig.l). Our 
explanation shows that the origin of the restriction comes from the boundary conditions. 



FIG. 1: The field patterns as a function of x and z for an ideal slab with e = = ~1. The arrows indicate the directions of 
energy flows. Inside the slab, energy flows in the direction opposite to that of the wave vector k. In (a) the slab has a thickness 
d shorter than the distance d\ — \zo\ between the source point and the left surface of the slab. The waves inside the slab region 
can be viewed of as radiated from a virtual source behind the slab, whereas the waves in the region behind the slab can be 
viewed of as radiated from a virtual source inside the slab. In (b) the slab has a thickness d larger than di = \zo\. No stationary 
solution can exist in the empty (question marks) region. 

Since the ideal slab does not reflect light at all Jj, thus the field inside and behind the the slab are 

Einsido{r,t)=AoH^o'\k\r-ri\)e-'^\ (13) 

and 

£;bchind(r,i) = Aoff^')(fc|r-r2|)e-*"*, (14) 

respectively. Here ri = (0, — zq) and r2 — (0, 2d+zo) — (0, 2d — di) are the positions of the two images predicted by the 
geometric optics. Now, if d < di, then ri and T2 are respectively located outside and inside of the slab, respectively; 
that is, they are virtual images (virtual line sources). In this case the fields are finite everywhere except at the source 
point. However, if we increase the slab thickness to d > di, then both images become real, and this contradicts the 
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boundary conditions. More specifically, a real image means a delta function term, i.e., a line source. Since there is no 
any other line source except the original one that located at rg, we conclude that the perfect imaging is impossible. 

Put it in another way. The time- averaged Poynting vector S must satisfy the divergenceless condition V • S = 0, 
thus there should be no singular point satisfying |S| = oo except the source point. Since in the slab the wave vector 
k and S are antiparallel to each other, thus the waves propagating in the < z < —zq and —zq < z < d regions must 
be "radiated from" and "absorpted by" the image inside the slab. This leads to the wave phase mismatch at z = — zq 
if ^0 7^ 0. From these consideration we conclude that the thickness limitation is a restriction originating from the 
boundary conditions of this system, and it implies that the stationary state (monochromatic waves) cannot satisfy 
these boundary conditions simultaneously. In other words, there is no stationary state. 

This result is consistent with the time domain results in Refill where the source was treated as a driving force and 
the two surface plasmon modes were two coupled oscillators. As one can see in Fig. 3 of Refill, the time evolution of 
the modulation amplitude A(t)/Astat{T — * oo) oscillates with a period Tosc = 47r/Act;fe, where Alu^ is the frequency 
difference between the symmetric and antisymmetric surface plasmon modes. When the absorption of the slab goes to 
zero, these two modes become degenerate, which leads to Tosc oo and AstatiT oo) — )■ oo. This case corresponds 
to the problem of driven oscillation without damping term. Therefore, the stationary state will not appear, and the 
field energy inside the slab grows to a larger and larger value without limitation. 

IV. A LOSSY SLAB 

Now we turn to the discription of the numerical results for a lossy slab. The permmitivity and permeability of 
the slab are chosen as e = —1 + iSe and ^ = — 1 + iS^j.; both Se and are small positive real numbers. With these 
parameters, the g function can be calculated. We then calculate the integral of Eq.Q) numerically as a sum. We 
first let kx = ktanO, with — /7r/2 < 9 < fn/2. Here is a reference angle, and < / < 1 gives the cutoff of 
{kx)max = (w/c) tan(/7r/2). In this paper we choose / = 0.96, which gives us a {kx)max/k « 16, large enough and 
numerically implementable to give us meaningful results about subwavelength imaging. The range (— /7r/2, /7r/2) is 
then being discrctized to = 3000 intervals, and the dk^ is replaced by ksec^{9)d9, with d9 — fir/ns- 

FIG. 2; (al) The field strength pattern as a function of x and z. In this case zq — —1, d — 2, \ — 0.3, e = /j. — —1 + O.OOli. 
The images have widths of the wavelength scale. The two straight lines represent the boundaries of the slab. (a2) The field 
strength at the focal plane as a function of x. (a3) The field strength on the x=0 plane as a function of z. The three straight 
lines represent the slab boundaries and the focal plane, (bl) to (b3) are for the case of subwavelength images. In this case 
20 = -1, d = 2, A = 2, and e = = -1 + O.OOOOOH. 

Figure 2. shows two typical cases for the imaging problem. In case A (Fig.2(al) to (a3)) the lens system creates 
two images, one inside and one outside of the slab, and they have widths of the wavelength scale. Here we have chosen 
zo = —1, d = 2, X = 211 /k — 0.3, and e = /i = — 1 + 0.001i. We observe clearly that the largest field strength locates at 
the two images. However, there is also some surface resonance effect near the boundaries. As we decrease the degree 
of the absorption, a stronger surface reresonance effect is observed. In case B (Fig.2(bl) to (b3)) we choose zq = —1, 
d = 2, A = 27r/fc — 2, and e ^ ji = — 1 + O.OOOOOli. In this case, the images become subwavelength scale. It is also clear 
that the field strength is very large at the two boundaries of the slab. This implies that surface-plasmon-polariton 
(SPP) plays important roles in this case. It is interesting to note that, although on the focal plane the field strength 
indeed has a peak along the cc-direction, however, the field strength does not have a local maximum around the image, 
and in the z-direction the wave strength decays from the second slab boundary. In this example the field strength 
at the focal plane is about only 1% of that at the boundaries. A closer observation find that the field strength at 
the image point is the same order as that around the source. This implies that if we turn on a line source, then the 
system has to spend a long time (several hundreds of 27r/w or above) to build the energy of the surface modes. Only 
after this transient process could the lens system focus the light to a subwavelength space region. 

The decaying profile of the field strength can be explained by the uncertainty principle. According to this principle, 
we must have the relation AzAfcx > 1, here the Aa; represents the width of the image, and the Afc^^ represents 
the fluctuation of fc^,. A subwavelength image is mainly formed by summing over the Fourier components of those 
\kx\ ^ w/c terms. Since k1 — uP- j (? — k\., these components must have imaginary /c^'s, and this leads to the decaying 
profile of the field strength. 

An approximate image size can be obtained by analyzing the transmission coefficient T. For \kx\ ^ w/c, and 
e = /x = — 1 + i^, (5<<1, we have 



(15) 
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which is a hyperbohc secant function with a peak value 



at the transverse wave number 

h = {lld)H2/5). (17) 

Thus the image size is given by 

A similar result has already been given by Merlin in Ref. l20l 

For the case B of Fig. 2, we have W ~ 0.43A, which is indeed a subwavelength focusing. However, the actual size of 
the image is in fact a little larger than that given by Ea. l|18|) . The reason is that the g{z^ kr^) function in the z > d 
region contains a factor e''"*^^"'*'^'^''! V(— 2iko), and thus in the integral (0)) the contributions from Fourier component 
with \kx\ < kx cannot be neglated. It seems that the near field excitations (evascent surface waves) and small enough 
absorption play the most important roles in the subwavelength imaging process. 



V. CONCLUSION 



In conclusion, we have studied the imaging properties of a negative-refraction slab lens, using a spectrum decompo- 
sition method. We have also given a simple geometrical explanation to the reason of the slab thickness limitation for 
an ideal negative refraction lens. For a slab with appropriate amount absorption, wc found that both the wavelength 
size and subwavelength size images can be formed. 
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